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Abstract 

We present a model of an economy inspired by individual based model ap- 
proaches in evolutionary ecology. We demonstrate that evolutionary dynam- 
ics in a space of companies interconnected through a correlated interaction 
matrix produces time dependencies of the total size of the economy total 
number of companies, companies age and capital distribution that compares 
well with statistics for USA. We discuss the relevance of our modeling frame- 
work to policy making. 
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1. Introduction 

Economies are highly complex systems. History suggests that the stan- 
dard analysis of an economy with the reductionist approach of individual ra- 
tionality and utility or profit maximization misses important features about 
aggregate dynamics and global stability that result from the interactions of 
economical agents [I]. Complexity theory offers an interesting alternative 
that has potential to provide better insights about the systemic risks by 
analyzing the properties of the economy as a whole [2] . 
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In this paper we develop and explore a simple model of an economy 
inspired by complexity models of evolutionary ecology, in particular the 
work by Laird and Jensen [3]|1]. This kind of analysis is relatively new 
in the economics literature, especially when it comes to policy design. In a 
world increasingly globalized with greater interdependence of the individual 
economies, analyzing systemic risks becomes critical in order to control and 
avoid global crisis. Andrew G Haldane, Executive Director of Financial Sta- 
bility department of the Bank of England, has published numerous papers 
about systemic risk of the financial sector in which he discusses interesting 
and innovative ideas for policy design Many of his arguments are 

based on complexity models of the global banking system where system pa- 
rameters, such as the density of interbank loans, have a drastic effect on the 
systemic risk. In general these complexity models consist of units interact- 
ing and forming networks. Important tendencies of the systemic risk can be 
drawn from properties of the network and the dynamics of the system such 
as the number of links or hierarchy of the network. Such properties could be 
taken in to account when designing economic policies and regulations. The 
first step to this aim is to establish a minimal modeling framework of the 
economy that allows us to reproduce, at least at a qualitative level, dynami- 
cal behavior consistent with real statistics, and that can be used to develop 
an appropriate concrete ecological approach to economics allowing us to rely 
less on purely biological analogies [7j[8j. 

We have modified the Tangled Nature model of evolutionary ecology de- 
veloped by one of the authors and his collaborators (see e.g. Jensen and 
Arcaute [9] and references there in). The reason we take our outset in the 
Tangled Nature framework is that it was demonstrated for the case of evolu- 
tionary ecology that this very minimalist approach compares very well with 
observation. From the simple assumption of mutation prone reproduction of 
individuals whose reproduction rate depends on the instantaneous configura- 
tion in type space, a long list of evolutionary and ecological phenomenology 
is reproduced. For example the intermittency of extinction and creation 
events termed Punctuated Equilibrium by Gouold and Eldredge (TU], realis- 
tic species abundance distributions and species area laws. The model can also 
reproduce relationships between strength of interaction and diversity seen in 
microbial experiments |TT]. This efficiency in obtaining systemic level phe- 
nomenology from some basic dynamical assumptions suggests to us that it is 
worthwhile to explore to what extend a similar framework can be developed 
for economics. 
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Since the foundation of the Santa Fe Institute almost three decades a go, 
researchers have explore the application of complexity science to approach 
economic problems avoiding the assumption that an economy is a system in 
equilibrium |ij . For a complete listing of research in complexity economics see 
Kling [12] and the references. Our contribution to this literature is to describe 
an evolutionary model of interacting components representing companies' 
interactions which aggregate into a model economy that agrees qualitatively 
with the evolution of the US economy since the great depression, a historical 
event in time defined by Kling [12] as a restructuring in the world economy 
that induced a new phase with a new economic structure, which we consider 
comparable with the beginning of our simulated economy. We find qualitative 
resemblance in the evolution of the Gross Domestic Product (GDP), the 
evolution of the number of companies in the economy, and the distribution 
of company ages and capital. In addition key parameters of the model that 
can be related to competition laws and the density of company business 
interactions and leverage, have important effects on the systemic risk of the 
simulated economy. We suggest that in the future our modeling approach 
may help to complement the standard equilibrium theory in macro economics 
by complexity economics. 

The remainder of this paper is organized as follows: Section 2 presents the 
model dynamics. In Section 3 we compare the behavior of our model to the 
US economy from 1929 to 2010. Section 4 presents preliminary predictions 
from the model regarding systemic risks of the economy with potential policy 
implications. Section 5 contains a brief discussion and concludes. 

2. The Model Framework 

We use a generalization of the version of the Tangled Nature framework 
described in Ref. [3j. A company a is represented by a string of L traits 
T" = (T{*, T^") in what we denote the economy space, and has capital 
C"(t) at time t. Here we have used G {0, 1, 2, 999} with periodic 
boundary conditions, i.e. there is only one unit distance between 999 and 
0. For computational and representational ease in this paper we have set 
L = 3. We can, for example, think of the three coordinates as indicative of 
the intensity in the use of inputs from the agricultural, industry, and services 
sectors which are embedded in final product of the company. Hence two 
companies close together in the economy space will produce similar products. 
Let us clarify that since the coordinates are periodical, in this interpretation 
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we need to consider the axis of coordinates as a circle and we can, say, think 
of the companies at the "top of the circle" (e.g. with ?=s 500) as producing 
a product highly specialized on the given sector. For example, an automobile 
company could be high on the industry coordinate circle, at medium hight 
on the services coordinate circle (e.g. it may have a financing department 
which lends to buyers, as well as publicity and others services involved), and 
low on the agriculture coordinate circle (e.g. they may use some wood in 
their cars). And for instance, say, a computing company will be high on 
the industry circle but it may be on the opposite side of the circle as the 
automobile company e.g. the International Space Station can be at the very 
top of the industry circle with coordinate 500, a car company can be at little 
lower at coordinate 550, and a computing company can be 450. There are 
1000'^ possible companies in the considered economy space. 

We use an stochastic sequential update for the capital C"(t) of company 
a at period i + 1 by following the update algorithm described below. One 
iteration corresponds to N{t) stochastic updates, where N{t) is the number 
of companies present at time t. The companies affect each other in two 
different ways: one is doing direct 'business', which we will call interactions, 
representing anything from buying or selling merchandise between them, to 
consulting and engaging in financial transactions; the other is competition 
and applies to companies that are close together in the economy space (i.e. 
produce similar products). The outcome of the interactions and competition 
together with the natural resources available in the system determines the 
probability of gaining or losing capital at the end of the iteration. 

Company a interacts with company /3 with strength J(q;, /3) independent 
of the corresponding reciprocal interaction strength for ^, J{f3,a). Only 
a subset of these a-p interactions are active (as in a real economy not all 
companies do businesses together) and the number of potential interactions 
depends on model parameter Cconnect (see Appendix A). We compute J {a, /3) 
by adding interaction contributions Ji{T°',Tf) from L different kind of in- 
teractions 



The individual interaction contributions are normally distributed so we use 
the jjj2 to ensure the total interaction is similarly distributed. For simplicity 
we use the same number of individual interaction contributions as the num- 
ber of coordinates, L = 3, and one can, for example, interpret these three 




(1) 



i=l 
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individual interactions as exchange of goods, financial transactions and intel- 
lectual collaboration respectively. In Appendix A we describe our procedure 
to obtain an exponentially correlated set of interactions J{a,P). 

An important property of these interactions is that two companies close 
together in the economy space have highly correlated interactions with any 
third company. By construction as the distance between a and P increases 
in the economy space there is an exponentially decay of the correlations 

Coruia,(3) = (jdTrX),MTf,T?')) (2) 



for any third company 7. This property has a natural interpretation. Say 
company a interacts with company 7 through trait i (e.g. exchange of goods), 
and say company P is very close to company a in the economy space (i.e. 
they produce very similar products), then Tf'^ will be very close to T"^ (see 



sum (A. 2) in Appendix A) hence mapping to coordinates close together on 
the corresponding trait axis, so the /3-7 interactions on trait i will be highly 
correlated to the a-7 interaction i.e. a and f3 will benefit (or be disadvan- 
taged) from company 7 with a similar strength provided that companies /3 
and 7 also do business on trait i (i.e. provided that the corresponding in- 
teraction series is active in the trait T^'^ so that the interaction Ji{T^ ^T^) 
is active). The degree of similarity of companies decays exponentially over a 
correlation length ^, a parameter of the model set at 300, which characterizes 



the matrix J(/3,a). For details see Eq. (A. 4) in Appendix A. 

Based on the property that the similarity of companies' profile decays 
exponentially with distance in the economy space, we define a direct compe- 
tition function C(a,/3) with a similar exponential decay, in particular with 
same decay constant ^ as the correlation length of the interaction series (see 
Appendix A for details). Companies close together in the economy space 
are similar due to the previous interaction correlation property, hence they 
compete strongly, whereas companies more than ^ units away in the economy 
space have negligible competition as their similarity also becomes negligible. 
We do acknowledge there may be indirect competition that arises from other 
economic sectors, that is from a company far a way in the economy space, 
which will not be capture by this competition function. This common in- 
direct competition are likely to be subject of technological innovation and 
evolution of the economy, and can be capture by the negative interactions 
J(/3,a). 

The performance of each company at each update is subject to natural re- 
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sources available in the system, R{t); each company takes one unit of natural 
resources from R{t) at its creation and when they go bankrupt they return 
their unit back to the pool of natural resources making the sum N{t) + R{t) 
a constant, where N{t) is the number of companies at time t. Obviously a 
more realistic description would be to let the resource bound in a company 
dependent on the type of company i.e. companies' position in the economy 
space. But for simplicity here we assume the resource to be independent of 
time and type of company. 

Using the total interaction strength, the total competition strength and 
the natural resources we can now define a weight function H{a,t) that will 
determine the probability of company a to gain or lose capital: 



where coefficients are just constants used to ensure the three terms have a 
similar order of magnitude. The first term represents the effect of interactions 
with all other companies present in the system at time t (note that many 
interactions are inactive, i.e. equal to zero). The denominator captures 
the number of companies in the same economic sector as company a. The 
rational for this ratio may be better explained by the following example: Say 
company a produces cars, then all other car companies 7 are points close to 
a in the economy space and due to the correlations of the interaction matrix 
the interactions will have similar total strength Yl^=i J{l^P)- But usually 
a given company P which is beneficial to company a can not realistically 
benefit all car companies (with some monopoly exception in the real world) 
hence to make it more realistic we divide the total interaction by the proxy of 
the number of companies in the given economic sector. The second term in 
the weight function ([s]) represents the direct competition affecting company 
a, therefore the term is negative capturing the negative impact of competing 
companies with the same profile and market aim. The third term represents 
the finite amount of resources available to the N{t) companies. The relative 
strength between these three terms is now used to determine the probability 
of gaining or losing capital. Namely, we define 





exp[H{a, t)] 



(4) 



gam 



1 + exp[i7(a,t)] 
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Upon success, company a's capital increases: 

C"(t + l) = C«(t)(l + c,-^) (5) 

where J'^{a) is the total strength of positive interactions, J'^°*{a) is the 
total sum of the absolute value of all of a's interactions, and Cg is simply a 
gain control coefficient to ensure realistic growth. Similarly, if unsuccessful, 
company a loses capital according to 

C'(t^l)=C'(t)(l-c.^). (6) 

A company whose capital goes below the bankruptcy threshold is simply 
removed from the economy and one unit of natural resources is returned to 
the pool R{t + 1). On the other hand if a company's capital goes above the 
investment threshold the company may invests in a new company with con- 
stant probability Pmv Upon success a new company is created in a random 
position in the region of size 100 x 100 x 100 centered around the position of 
the founder company in the economy space. The new company starts with 
10% of the founder's capital and takes a unit of natural resource from the 
pool R(t + 1). Note that the new company will be well within correlation 
length ^ (fixed at 300) form the founder company. This means that the new 
company will have a negative impact on its founder in terms of competition - 
although there can be a trade-off in terms of positive interactions. Therefore 
this feature of the model is to be interpreted as a competition law within a 
successful market rather than an expansion of the investing company. For 
example think of when the US supreme court divided Standard Oil in 1911 
because they were dominating the oil market. We will see that the proba- 
bility of investment Pinv is hence a very important parameter as it controls 
the distribution of capital and the appearance of oligopolies in the system. 
We will see in Section 4 details of the impact on the systematic risk of the 
economy. 

We let the model run for several thousand iterations and observe the 
aggregate dynamics that arise from these simple companies interactions. 

3. Case Study 

To assess the usefulness of our simple abstract approach we now make a 
comparison between dynamical trends in the model and equivalent observ- 
ables from the the US economy, where data can be readily obtained. The 
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model is able to generate aggregate dynamics that are similar to those ob- 
served in the US economy. The dynamics of macroeconomic variables are 
related to key model parameters such as the number of potential business in- 
teraction in the economy space ( Cconmct), and the redistribution of capital by 
the probability of generating new companies {Pinv)- Nonetheless contingency 
plays a very important role: two independent runs with the same parameters 
typically exhibit evolutionary paths that differ in details; but this is certainly 
consistent with experience from the real world. 

General Macroeconomic Trends: Model v.s. US Economy 

Our starting configuration consists of 1000 companies placed at random 
positions in the economy space and with a randomly assigned capital in the 
range between 110% of the Bankruptcy Threshold and 110% of the Invest- 
ment Threshold. After about 1000 iterations this random set of companies, 
with a given matrix of interactions, self-organizes into well defined structures 
in the economy space consisting of tubes and planes (see Figure ([T|). These 
structures can be interpreted as the supply chains of given products. Like 
in real economies who specialize in specific productions, the simulated econ- 
omy also undergo specialization, which is seen as the formation of coherent 
structures around specific trait coordinates. Hence, probably not controver- 
sial at all, our model indicates that specialization of economies is mainly an 
emerging property driven by increasing returns and path dependence. 

Here we focus in the following macroeconomic variables among other fea- 
sible observable variables: GDP (as the proxy for GDP in our model we refer 
to the sum of the capital of all companies at time t), the total number of 
companies, and cross section distributions of companies' age and capital. We 
ran 40 independent simulations with the same model parameters but random 
initial companies, out of which 15 simulations collapsed after a few iterations 
- they did not survive the initial chaos. We are interested in the average 
behavior of the remaining 25 simulations. All graphs referring to the model 
represent averages of these 25 simulations. The GDP growth rates are com- 
puted with a lag of 100 iterations giving realistic values for annual growth 
rates. We use this comparison to calibrate the time scale of the model against 
real economical time i.e. think of 100 iterations of the model as represent- 
ing one-year period in the real world. We use US data from 1929 to 2010; 
this period of time for the US economy provides a good comparison since it 
starts at the beginning of the great depression which brought the economy 
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into chaos and marks a 'new beginning', a new restructure much hke the 
start of the simulated economy. 

We find that both the US and the simulated economy have a volatile 
economic growth at the beginning of the period, with periods of fast GDP 
growth followed by recession. As both economies mature we observe a damp- 
ing evolution on their GDP growth rates towards a more modest and stable 
growth (see Figure ^ and (|3])). The evolution of the GDP is similar for both 
economies, both enjoying an overall growth (see Figure IQ and (g), but an 
important difference is that the US GDP follows a convex evolution whereas 
the simulated GDP tends to be concave. The reason for the concave shape 
of the simulated GDP comes from the effect of the fixed natural resources 
limit which restricts the economic growth. 

Another interesting comparison is the distribution of companies size in 
the model and in the US economy. Figure ^ shows a histogram of US firms 
in 2007 sorted by receipts size (i.e. expenses) together with model results 
for the companies capital distribution across the economy at t = 49 years 
(note this is already an steady period)]^ Again we notice nice similarities, the 
model too has many small companies and a few big ones and the proportions 
are realistic. 

We also compare companies' age distribution across the economy. Figure 
([7]) compares histograms of the US firms by age in 2007 to average age dis- 
tribution at t = 49 years. They both have a decreasing number of companies 
as the age categories increases. Note that we ignore the first age category in 
the model (less than 75 iterations old) because most of these new companies 
in the model only live for a few iterations. We argue that these companies 
would not make it to the 'national statistics'. But indeed one of the limita- 
tions of our model is that it has unrealistically large creation and destruction 
rates (some times reaching as much as 200% per 'year' !). This property can 
be tuned by manipulating the thresholds for bankruptcies and for creation 
of new companies through investments. 

Finally the evolution of the number of US firms from 1977 to 2010 and 
the evolution of the number of companies in the model are shown in Figure 
(jsj) and (|9|. The data of the US starts in 1977 hence for consistence with 



^We chose 2007, and accordingly t=4900 in the model, in order to avoid potential 
unrepresentative perturbation in the US data by the finantial crisis that started in 2008. 
Nevertheless there is no qualitative difference when comparing 2009 to t=5000. 
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previous interpretations we compare with the last 2500 iterations of the model 
where we have a more stable growth similar to the one observed in the US 
data. The model has a similar tendency although it has more variance due 
to the problem of excessively large companies destruction and creation rates 
previously mentioned. 

4. Systemic risk of the model 

We would like to explore the potential use of this model. To this end we 
analyze how the model dynamics change when we modify two key parameters: 
The first one is the probability of investment Pmv, which is related with the 
distribution of capital in the economy, much like a fair competition law; 
the second parameter is the connectivity threshold of the interaction series, 
C connect, which is related to the number of potential interactions between 
companies in the economy space and may be related to a confident / reheated 
economy for high connectivity, and vice-versa, and also related to companies 
leverage regulations. We run several simulations for the following values of 
Pinv = {0.1; 0.3; 0.5} and with Cconnect = { — 10; 0; 10} (note that in this part 
of the analysis we did not use averages of multiple simulations due to time 
constraints, but the patterns we will discuss are very consistent within the 
simulations we ran for each parameter combination, that is at least three 
simulations per parameter combination)]^ To quantify stability and risk of 
the system we consider the mean and the volatility (i.e. the variance) of 
the growth rates of macroeconomic variables: GDP, number of companies, 
average capital per company, and companies creation/destruction. Due to 
space constraints we will not go into details of the results, but only discuss 
the volatility of GDP to describe the general behavior of the model. 

We found that an economy with loose fair competition regulations (i.e. 
low probability of investment Pinv) increases the average GDP growth but 
also increases the volatility of this variable inducing strong recessions in the 
simulated economy. Whereas high probability of investment generates a more 
modest and stable GDP growth. Similarly a reheated economy with loose 
leverage regulations (i.e. highly connected economy) tends to have a strong 



•^Recall from the model dynamics in Section 1 that lower connectivity threshold gives 
rise to high number of potential interactions in the system. Hence we will refer to scenarios 
with Cconnect — ~10 as high connectivity, Cconnect = as medium connectivity and 
Cconnect — 10 as low Connectivity. 
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and volatile GDP growth, and low connectivity induces modest but stable 
growth. As we test different combinations of these two parameters we find 
what we would expect: deregulation in competition laws in a reheated econ- 
omy i.e. low probability of investment in a highly connected economy, is the 
most unstable combination; and strong competition regulation together with 
leverage regulation i.e. high probability of investment in an economy with 
low connectivity, is the most stable and modest growth. Further more, the 
unstable parameter combinations tend to develop oligopolies in the economy, 
an small number of companies that completely dominate given economic sec- 



tors (see Figure ( 10 ) ). This explains partly why these scenarios are unstable, 
as the relevant economies are highly dependent on few companies, whereas in 
the economies which favor a fair competition, that is the stable scenarios, the 
economic growth relays on many companies so that the collapse of a single 
company is almost unnoticed. 

This results are not surprising. With low probability for investment com- 
panies can grow much bigger, and some eventually become oligopolies or even 
monopolies with no regulation to control them. This is further facilitated in 
a highly connected economy were key positions will be at the center of a 
large web of business, hence companies at advantageous locations are likely 
to become very successful, although still liable to chance, and the fate of 
the economy ends up in the hands of few companies. High conectivity also 
makes the economy more interdependent and the bad fate of few companies 
can expand trough the economy with easy generating economic recessions. 
Hence arguably we could identify the combination of unstable parameters 
as corresponding to a deregulated economical policy. Both in terms of lack- 
ing fair competition regulations (low probability of investment), and allowing 
high leverage for companies and deregulated banking system in terms of high 
number of loans and debts (high connectivity)]^ Our model is still to simplis- 
tic to provide policy advice but is very encouraging that it agrees with well 
known risk of deregulated economies and furthermore it captures very well 
the tradeoff between higher average growth and higher systemic risk from 
deregulation. 

In the future such models could be tuned with real data from an economy 
taking in account major agents such as big banks and corporations in order 



"^Note that our result concerning the effect of connectivity on the systemic risk of the 
economy is consistent with the finding of A.G.Haldane and May [6]. 
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to start the simulation with the present state of the economy. Prom there we 

could run several simulations and test several policy strategies, and quantify 
the performance of potential policies in a way that classic economics, in 
particular macro economics and the general equilibrium theory, cannot. 

Discussion 

The simple model developed here is able to generate aggregate dynam- 
ics for key economic variables that are generally consistent with empirical 
evidence. There are nonetheless several limitations and opportunities for im- 
provement. One issue is the effect of the natural resources limit which tends 
to give a concave shape to the evolution of the GDP and of the number of 
companies in the model, in contrast with the convex shape shown by the US 
GDP. We notice that the natural resource limitation of the US is less of a 
barrier since on top of their national natural resources they can also import 
from other countries. We could design extensions of the model modifying the 
role of natural resources in the dynamics to try to capture better this fact 
and to allow for technological innovation. 

As in the real world contingency plays also an important role in the 
model: two independent runs with the same parameters typically exhibit 
evolutionary paths that differ in details. This emphasizes the importance of 
having adaptive policies that are able to introduce real time corrections when 
a given economic system starts to deviate from a stable path. The other im- 
portant observation in the model relates to patterns of specialization. Stable 
economies are not good at producing all goods. Gompanies end-up operat- 
ing in given regions of the economy space suggesting that the comparative 
advantages of an economy are not an exogenous but an emerging property 
in an economy. Finally, the model predicts the emergence of oligopolies and 
monopolies as a natural feature of deregulated economies. This also calls for 
policy interventions that regulate the level of competition, such as anti-trust 
law. 

We believe the importance of our model is that it highlights how the 
macro behavior, i.e. the growth or decline of companies, may originate in the 
entangled nature of the interaction between companies. We have deliberately 
kept the internal structure of the individual companies very simple in order 
to focus on the fact that a given company always finds itself embedded in an 
ever changing environment composed of all the other co-existing companies. 
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This fact suggests that the notion that a company may possess an intrinsic 
robustness or fitness is to an extend an illusion. 

In the future we plan to study the dynamics presented here on econom- 
ical networks extracted from real data sets of the network of transactions 
between companies. With increasing computer power and more sophisti- 
cated complexity models these approaches may become key in the future of 
macroeconomic policy design. 
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Fig. 1: Economy space at t=5000. Companies size proportional to age. 
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Fig. 4: US GDP 1929-2010. GDP corrected for inflation. Source: Bureau of Economic Analysis. 
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standard deviation divided by 10 away from the averaged given by the solid curve. The fluctuations are 
signiflcant. *Wc use 1 Year = 100 Iterations, see the text. 
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Appendix A. 

To compute the individual interactions Ji(T°',Tf) we use a method de- 
veloped by Laird and Jensen [3]. We use L independent autoregressive time 
series, each one associated to an individual interaction: 



_1 = pXn + 



(A.l) 



where is a normal distributed random shock and p is a constant that 
determines the correlation length of the series. We then map the following 
sum to the corresponding time series 



(A.2) 



k=l 



where bij,Cij G { — 1,1} are independent randomly assigned and kept fixed 
trough out the dynamics. The individual interaction i will be the value of 
the time series corresponding to the index T^'^ (i.e. Ji{T°',Tf) = Xr^ap be- 



longing to the i series). We will refer to these three time series as the set 
of interaction series. To create the subset of inactive interactions, for each 
interaction series, we create a parallel independent 'switch' time series and 
whenever this second time series goes below a given connectivity threshold, 
Cconnect, the Corresponding values of the interaction time series are set to zero 
(see Figure ( A.ll[ )). This connectivity threshold is important as it will influ- 
ences the number of potential interactions in the economy (i.e the number of 
connected points in the economy space). 



The motivation to use sum (A.2) and the auto-regressive time series for 



the interactions comes from the properties of the autocorrelation of the time 



series: From equation (A.l) using induction we can write 



pXn-l + 

p{px 



n-2 



+ 



+ 



P Xn-2 + P0n-1 + 



Xr. 



fe=0 
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Fig. A.ll: Interaction series (top) and corresponding switch series (bottom) with C connect = 0. 

assuming xn = and taking — )■ — oo. Now are IID normal N(0,1), 
therefore clearly we have E{xn) = 0, also cov{(f)i, (pj) = for i 7^ j and 
E{4>'^_^) = var{(j)n-k) = 1 hence we have: 

00 ^ 
var{xn) = E{xl) = p"'E{<Pl_,) = — ^ 

fc=0 ^ 

provided p < 1 (which we ensure in the model). Also p < 1 implies that the 
process is stationary hence the auto-covariance and autocorrelation between 
Xn and Xm depend only on the separation length r = \n — m\. Define = 
cov{xn,Xn-T), we have: 

Xn = PXn-1 + 0n 

taking expectation 

S-f E{XnXn—T) pEiXn—lXn—r) 
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= PSr-l 



P So 



using induction and the fact that Xn-r depends only on (pn-r and earher 
shocks therefore Xn-r is uncorrelated with latter 0j's (i.e. E{(j)nXn-T) = 
for r > 1) . Finally the most important property of the AR(1) process for 
the purpose of this model is the autocorrelation sequence which takes the 
following form 

Corri^Xn, Xn-r 



So 



p 



Hence the autocorrelation of the time series has an exponential decay 



Corr{xn, Xr 



\n—m\ 



P 

exp 



(A.3) 



\n — m\ 



For the purpose of this model we choose the correlation length ^ and we set 

p = exp[-l/^] 

Since the interaction behavior will be considerably correlated for compa- 
nies closer than ^ units away on each coordinate, we consider these companies 
as producing similar products, hence the correlation length represents the 
competition length in the economy space i.e. the distance between compa- 
nies at which they should have a non-negligible competition between them. 
But to define the competition strength, instead of using correlations for each 
interaction series, we define a single function representing an exponentially 
decaying competition strength that approximates a similar two point corre- 
lation function determined by the separation of the two companies in the 
economy space: 



C{a, (3) = exp 



-l/L)Eti AT, 



a/3 



(A.4) 



where 



AT;"^ = |6,i(Tf - Tf) + k2{T^ - T^) + ... + hLiTt - T, 



and hij are the same we use in sum ( A.2 ) . Self-interactions and self-competition 
are dismissed. 
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